


X-551-72-277 


TROPOSPHERIC RANGE-RATE 
TRACKING-DATA CORRECTION 


John W. Marini 


August 1972 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 

f 

I 



TROPOSPHERIC RANGE-RATE 
TRACKING- DATA CORRECTION 

John W. Marini 


ABSTRACT 


A formula for correcting the tropospheric error in range- rate 
satellite tracking data is given. The formula is based on the 
method, already published, which was used to obtain correc- 
tions for elevation-angle and range data. 

In addition, an improved method is given for calculating some 
of the parameters required in the correction formulas. 
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TROPOSPHERIC RANGE-RATE 
TRACKING-DATA CORRECTION 


INTRODUCTION 

The corrections given here are based on a method already published [1, 2] . 

The method employs a ratio of polynomials that is adjusted to agree simultane- 
ously with two different series expansions, one of which holds at small values 
of the elevation angle, while the other, an asymptotic expansion, holds at large 
values of the angle. 

Most satellite tracking is done at elevation angles above about 10 degrees. At 
these higher angles, the magnitude of the tropospheric correction depends 
mainly on the coefficients of the terms in the asymptotic expansion, and these 
may be calculated or estimated provided only that the surface refractivity and 
the zenith integral of the refractivity from the ground station up is known. The 
surface refractivity can be calculated and the zenith integral can be estimated 
from measurements of pressure temperature, and relative humidity, taken at 
ground level at the tracking station. 


CALCULATION OF SURFACE REFRACTIVITY 

The surface refractivity N 0 can be calculated from the surface pressure P 0 
(millibars), the temperature T (degrees Celsius), and the relative humidity R h 
(percent). The formula for refractivity is [3] 


77.6P 3.73 x IQ 5 e 

T + 273. 15 + ( T + 273.15) 2 


where e (millibars) is the partial pressure of the water vapor, which can be 
calculated from the relative humidity and temperature using [ 4] 


e 


R h 

100 


7.5 T 

X 6.11 X 10 237 - 3 + T 


(2) 


Alternatively if the wet bulb temperature T w (°C) of a psychrometer is given 
instead of the relative humidity, e may be calculated from [3,5] 


1 



7.5 T 


(3) 


w 

e = 6.11 x 10237.3 +TF^ _ 0.00067 P (T - T w ) 


Finally, if the dew-point temperature T d (°C) is given 


7.5 T d 

e = 6.11 X 10237.3 + Td 


(4) 


ESTIMATION OF THE ZENITH INTEGRAL 
The zenith integral is defined as 



where N (h) is the refractivity at the height h (km) above sea level, and h 0 is 
the height of the surface where the tracking station is located. 

This integral is evaluated by separating the refractivity N into two parts 


N - N d + N w (6) 

where N d is the contribution to N from the polarizability (nominally the dry 
part) of the air and N w is the contribution from the (dipole moment of the) water 
vapor that is present. The zenith integral becomes 


^ zd ^ ^ zw 


|-*00 s*oo 

N d dh + N w 

•'hn J hn 


dh 


(7) 


The dry part may be evaluated from a knowledge of surface pressure P Q (mb) 
[6,7] 


I 2d = kP 0 (km) (8) 

where k equals 2.2757 at 45° latitude and varies with latitude as shown in 
Figure 1. 
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PARAMETER k 


Reproduced from Reference 7 


£ OBSERVED ± a (a OMITTED FOR VERY CLOSE POINTS) 

A CALCULATED USING g AT OBSERVED ATMOSPHERE CENTROID, 
APL GEODESY, DRY AIR 


CALCULATED FOR ELLIPSOIDAL EARTH, DIFFERENT CENTROIDS, 
DRY AIR 



LATITUDE (degrees) 

Figure 1. Latitude Variation of Tropospheric Range Error Parameter k 
(Where /N d dh = k P s ) 
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The wet part is not so easily evaluated, since there is no simple theoretical 
relationship between the integral of N w and the surface pressure of the water 
vapor in the air. Fortunately the contribution of the wet part is somewhat 
smaller than that of the dry, so that the relative error in N caused by inaccura- 
cies in N w is not as serious as it might otherwise be. As an interim measure, 
the wet part of the zenith integral may be estimated from 


N 


Ow 


f 


(N w /N 0w )dh - N 0w H w 


2 N 0 w ( km ) ( 9 ) 


where N 0w is the surface value of the wet part of the refractivity (the second 
term in Eq. 1), and the factor 2 arises from assuming an equivalent height 
(h w - h 0 ) [6] of 10 km, or an exponential scale height H w of 2 km for the wet 
term. Work is in progress to provide a better estimate of I hw at each tracking 
station through statistical analysis of appropriate radiosonde data. 


CALCULATION OF THE SCALE HEIGHT 

The correction equations do not employ the zenith integral I directly, but are 
formulated instead using the scale height H defined as 


H = I z /N 0 (10) 

with the surface refractivity calculated from (1) using (2), (3) or (4), and the 
scale height calculated from (7), (8), (9) and (10), the calculation of the correc- 
tions is a straight-forward application of the equations already published [ 2] * 
together with the range- rate correction derived in the next section. 

The programs in APL language [ 8] are given in Appendix 1. 


DERIVATION OF RANGE-RATE CORRECTION 

The range- rate correction AR can be obtained by differentiation of the range- 
correction. If (43) of reference [2] is differentiated with respect to time, and 
if the dependence of the correction on R is neglected, there results 


The method for calculating one of the polynomial ratios has been modified. This change is 
described in Appendix 2. 
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A R 1 

— — = t: 10 -6 N Q cos d 0 [M'(a) - p QL(a) L'(a)] 6L (11) 

r o z 

Here M(a) has the expansions 

M(a) - (1/a) - Mj ( 1/ a) 3 + M 2 (l/a) 5 - . . . (12) 

M(a) = m Q - mja + m 2 a 2 - . . . (13) 

where the values of Mj, M 2 , m 0 , and m 1 are given in [2] and m 2 has the value 

m 2 = io ^ + i oO + q ^1 + | qi 0 2 j i 2 (14) 

and where 

i 2 - \ I"(0) 

has, for an exponential profile, the approximation 

i 2 = Vtt (1 - 1.023 q)" 1 - 8 ± 0. 15% 0 < q < 0.7 (15) 

Consequently M' (a) in (11) has the expansions, obtained by differentiating (12) 
and (13), 

NT (a) - - (1/a) 2 + 3Mj(l/a) 4 - 5M 2 (l/a) 5 + ... (16) 

M'(a) = - m 1 + 2m 2 a - . . . (17) 

which can be modelled as shown in Appendix 6 of [ 2] . Greater accuracy, how- 

ever, is obtained by using the procedure given in Appendix 2 of this report. 

The faction L'(a) in (11) equals 
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L'(a) - - a I'(a) - 1(a) + - q 1(a) I'(a) 


in which I' (a) can be modelled by the same procedure used for M'(a). 
Similarly, (67) of [2] gives 


AR l 


\ 10-6 N o E W - Qp (l - /3U - \ qu 2 ) • (u + J3\J' + \ qUU') E 


and W ' can be modelled in the same way using the result that 


J 2 ~ 2 w "(°) = ( u o 1 ~ \ q u i) 


Similarly, U'(/3) is modelled using Equation (65) of [2] which should read 


u 2 = ? U"(0) 


\ r (| q u o)/ 1 - \ qi’ (| qu 0 
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APPENDIX l 

APL PROGRAMS FOR CORRECTIONS 

Inputs to PREPASS 1 

Surface refractivity NNO 
Scale height HH (km) 

Distance RO (km) of tracking station from the center of the earth 
Outputs of PREPASS 1 

The 16 constants J1 - 14, M 1 - M4, SI - S4, and IP1 “ 

Inputs to CORR 1 

The 16 constants 11 - 14, Ml - M4, SI - S4, and IP1 - IP4 from PREPASS 1 
Surface refractivity NNO 
Scale height HH (km) 

Arrival angle THETAO (mrad) of the wave from the satellite. 

Arrival angle rate THETAO DOT (mrad/ sec) 

Range RR (km) of the satellite from the tracking station. 

Outputs of CORR 1 

Elevation angle correction DEE (mrad) 

E .= 9 0 - DEE 

Range correction DRR (km) 

E = - DRR 

Preceding page blank 
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Range-rate correction DRRDOT (m/s) 


R 


O 


R 


m e 


a s u r e d 


- DRRDOT 


Inputs to PREPASS 2 

Surface refractivity NNO 
Scale height HH (km) 

Distance RO (km) of tracking station from the center of the earth. 

Outputs of PREPASS 2 

The 16 constants U1 - U4, UP1 - UP4, W1 - W4, and II - 14. 

Inputs to CORR 2 

The 16 constants U1 - U4, UP1 - UP4, W1 - W4, and T1 - T4 from 
PREPASS 2. 

Surface refractivity NNO 
Scale height HH (km) 

Elevation angle E (milliradians) of the satellite above the tracking station 
Elevation angle rate EDOT (mrad/sec) 

Range RR (km) to satellite from the tracking station 

Outputs of CORR 2 


Elevation angle correction DEE (mrad) 



Range correction DRR (km) 


R 


R 


m e 


a s u red 


Range- Rate correction DRRDOT (m/ s) 

o ° 

R = R 

measured 


- DRR 


DRRDOT 
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VPREPASS1 [ □] V 

V PREPASS 1 

Cl] PSQUARED+-2*HHiR0 

[2] P+PSQUARED* 0.5 

[3] Q4-lE-6*NN0*R0iHH 

[4] F0«-I0«-( (Ol)*0.5)x(l-0.9206xQ)*"0.4468 

[5] F1+-I1+2 *1-Q 

[6] FFl-«-JJl-*-0 . 5*1-0. 5x<2 

[7] FF2<-I I2+-Q .75xl+$x“0.75+$t6 
C 8 ] CALCULATEF 

[9] Il<-PSQUARED*ll 

[10] 12+PSQUARED*F2 

[11] l3-*-PSQUARED*F3 

[12] J4+-PXF4 

[13] K 0+( (02 )*0„5)x(l-0.9408xQ)*" 0.4759 

[14] F0-*-M0+- ( J0xl+Qxi+gxJ0xJ0*12 )-G x F0t2 

[15] Fl-«-Ml-«-2x ( 1 + QxJ0xJ0*4)t1-<2 

[16] FFI+MM1+0 . 5x1-0 . 75x$ 

[17] FF2<-A/A/2^0.75xi+ex( _ 2 5f24) + llxg^36 

[18] CALCULATEF 

[19] Ml+PSQUAREDxFl 

[20] tt2*-PSQUAREDxF2 

[2 1] M3*-PSQUARED*F3 
[22] ^4^P«F4 

[2 3] J2-*-((oi)*0. 5 )x( l-i .02 3*fl)*"l . 8 

[24] M2-*-(I0 * ( 1 -Q ) *2 )+fixJ2xi+o .25x£xJ0*2 

[25] G0+M1 

[26] G\+-2xM2 

[27] GGl+ZxMUl 

[28] GG2+S*MM2 

[29] PREPAREF 

[30] CALCULATEF 

[31] £1*-PSQUARED*F1 

[32] Z2+PSQUARED*Z2 

[33] S.3*-PSQUARED«F3 

[34] S4«-Px£4 

[35] G0+I1 

[36] G1+2*I2 

[37] GGl+3'xII 1 

[38] GG2+5*II2 

[39] PREPAREF 
[4 0] CALCULATEF 

[41] IP1+PSQUARED*F1 

[42] IP2+PSQUARED*Z2 

[43] lP3*-PSQUAREDxZ3 

[44] IP4«-Px£4 

V 


12 



vcoj?j?i[[]]v 

V CORRhSIN iCOS‘,1 1 LL 1 M 1 S 

[1] SIN-loTHETA0*0 . 001 

[2] C0S—2OTHETA 0x0.001 

[3] I-iSIN+Il *SIN+12 iSIN+Z 3 *SItf+£4 

[4] LL+l-IxSIN-O . 5xlE~&xriN0*I 

[5] DEE-0 .OOl*NNOxCOSxI-ROxi,L*RR 

[ 6 ] M-*SIN+Ml*SIN+M2iSIN+M3iSIN+MH 

[7] DR R - 1E~&*NN0 xRH *M-0 . 5 x lF - 6 xAWO x ( (R0*C0S*LL) *2 ) *RR*HH 

[8] S-iSIN+SliSIN+S2iSIN+S3*SIN+£k 

[9] IP-( * SIN+IP1 *SIN+IP2 iSIN+£P3 tSIN+IP 1 * ) *2 

[10] TERM-Q*(R0iRR)x(.COS*2 ) *LL* J + ( 1F _ 6 *NN0* lx IP ) -SIN* IP 

[11] DRRDOT-~lE~9*THETAODOTxNNOxHHxCOSx(,S*2 ) - TERM 

V 


VPREPAREFlUlV 

V PREPAREF 

[1] F0-G0*0.5 

[2] Fl—Gl *2 *F0 

[3] FF1-GGH2 

[4] FF2-0 . 3*GG2 -FF1*2 

V 


VCALCULATEFWV 

V CALCULATEF 

[1] Fl—FFl 

[2] l2-(FF2iFl)-Fl 

[3] F3-*-F2t(F0xF0xF1x1+F1tF2)-1+F1xF1 

[4] F4«-F0x£1xF3tF2 

V 
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VPi?£P,4SS2[[]]V 
7 PREPASS2 

Cl] PSQU ARED*2 *HH *R0 

[2] P*-PSQUARED*Q . 5 

[3] Q4-1E~&*NN0xR0*HH 

[4] IPU-*-~2* ( 1 + 1 .482x0*” 0. 3826 
C5] DEN*- 1-0 . SxQxJPU 

[ 6 ] IPPU*-2x ((oi)*0.5)x(l + 1.71xfi)*0.1 

[7] F0*UQ*( (Ol )*0 .5)x(i+i.4844*Q)*“0. 39144 

[ 8 ] F1*-U1*--IPU+DEN 

C9] FF1*-UU1*0 . 5xi+0 . 5*Q 

[10] FF2*-UU2*-0 .75xl+$x(7*i2 ) + (J * 6 

[11] CALCULATEF 

[12] gl*PSQUARED*Fl 

[13] U2*-PSQUAREDxF2 

[14] U3*-PSQUAREDxFZ 

[15] £/4*-Px£4 

[16] G0*U1 

[17] Gl*2xU2*-0 . S xippu *DEN* 3 

[18] GGl*-3xUUl 

[19] GG2*-$xUU2 

[20] PREPAREF 
[•21] CALCULATEF 
[22] UPU-PSQUAREDxFl 
[2 3] UP2*-PSQUAREDxF2 
[24] UP3*-PSQUAREDxFZ 
[2 5] i/P4«-PxF4 

[26] KU*-( ( 02 )*0.5)x(l + l.6454xQ) *“o.583 

[2 7 ] FO*-W 0 "*- ( UOxl+Qxi-QxU0xUO*6)-O.5xQxKU 

[28] Fl*-Wl-*2 xl-0 .25xQxU0x[/0 

[29] FFl*-WWl-*0 .5x1 + 0. 25 X Q 

[30] FF2*-WW2+0 . 7 5 x l+Qx ( 7 * 24 ) +Q i 16 

[31] CALCULATEF 

[32] W1*-PSQU AREDx Fl 

[33] W2*-PSQUAREDxl2 

[34] W3*-PSQUAREDxl3 

[35] f74«-Px£4 

[36] J/2+i/0xl-0 . SxQxUl 

[37] G0*-W1 

[38] Gl-*2xW2 

[39] GGl*-3xWW 1 

[40] GG2*-5xWW2 

[41] PREPAREF 

[42] CALCULATEF 

[43] Z1+PSQUAREDxF1 

[44] Z2*-PSQUARED*Z2 

[45] ZZ+PSQUAREDxFz 

[46] £4«-Px£4 
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VCORR 2[G] V 

V CORR 2 iSINEiCOSE iU-.SRUPiPAREN iW-.T 

[ 1 ] SINE+IOEE* 0.001 

[ 2 ] C0SE+2OEEX0 .001 

[ 3 ] U<-*SINE+U1*SINE+I12*SINE+!13*SINE+UH 

[4] £/P«-( iSINE + UPliSINE + UP2iSINE + UP3iSINE+UPii )*2 

[5] PylPFAN-l-PxSItfp + O . 5xlp“6x//P0xj/ 

[6] DEE+-Q . 0 01*NN0x CO SE*U -PARER *( 1 - 1P _ 6 xtfffo x£/p ) xPO *PP 

[7] kV^5JArp+ViiSI//P + i/2T5J//p + ^3+SIWP+V4 

[8] DPP-*- IP - 6 *NNO xHH*W+Q . 5 x l£ _ 6 xtftfo x ( t H H*RR ) * (RQxCOS E*PAREN ) *2 

[9] T+-iSINE + £HSINE + T2i5INE + T3*SINE + Z^ 

[10] TERM+Q* £P 0 *PP )x ( COSE *2 )xPARENxU- ( 1P~6 xRROxUxUP)+SINExUP 

[11] DRRDOT+ IE~9xEDOTxNNOxHH*cOSEx(T*2)+TERM 

V 


15 



PRECEDING PAGE BLANK NOT FILMED 

APPENDIX 2 

APPROXIMATION TO G(a) 


Let the function G(a) have the expansions 

G(a) - (1/a) 2 - Gj (1/ a) 4 + G 2 (l/a) 6 ... (2-1) 


G(a) = g Q - g x a + . . . 


(2-2) 


for a large and small respectively. The modelling of G(a) is accomplished by 
approximating its square root as a ratio of polynomials. The square root of 
G (a) has the expansions 


VW) - (1/a) - \ Gj ( 1/a) 3 + \ G?j 


(1/a) 


5 - 


(2-3) 


/G^) 



1 !l 

2 fZo 


(2-4) 


as may be verified by squaring (2-3) and (2-4) to obtain (2-1) and (2-2). Expan- 
sions (2-3) and (2-4), however have the form of (38) and (39) of reference [2] , 
and so may be modelled by (33) of [ 2] . 

Hence 


G(a) = F 2 


2 G i 


G 0 ” „ 


Gi 2 


VT 0 1 gi / 2 


where F is defined by (33)-(37) of [2], 
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